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ABSTRACT 

It  is  shown  that  an  adaptive  multi-input,  multi-output  space-time 
processor  (beamformer  and  filter)  consists  of  a  multi-input,  multi¬ 
output  adaptive  processor  and  a  nonadaptive  beamsteering  operator. 

The  adaptive  processor  is  determined  by  the  spatial  and  temporal  prop¬ 
erties  of  the  noise  field;  whereas,  the  beamsteering  operator  is 
determined  a  priori  by  the  desired  beamsteering  directions.  In  this 
structure,  the  multichannel  adaptive  processor  is  common  to  all  the 
desired  beamsteering  directions. 
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I.  INTRODUCTION 

A  number  of  papers  on  the  adaptive  beamforming  literature  have 
considered  a  multi- input,  single  output  adaptive  beamformer  as  shown  in 
Figure  1  [l]-[6].  The  comment  is  usually  made  that  this  postulated 
structure  can  be  used  to  study  adaptive  beamformers  that  have  off- 
broadside  beamsteering  directions  by  placing  a  beamsteerer  between  the 
sensors  and  adaptive  beamformer.  The  difficulty  with  this  approach  is 
that  for  each  new  beamsteering  direction,  the  adaptive  filter  coefficients 
must  be  recomputed.  If  multiple  simultaneous  beamsteering  directions 
are  desired,  simultaneous  adaptation  for  all  the  beamsteering  directions 
becomes  intractable. 

An  alternative  approach  to  the  previously  cited  work  is  not  to 
postulate  a  structure  in  advance,  but  to  let  the  space-time  processor 
structure  evolve  from  the  mathematics  for  the  optimum  processor.  This 
approach  has  been  taken  by  Stocklin  [7],  Bryn  [8],  Middleton  and 
Groginsky  [9],  Adams  and  Nolte  [10],  Cox  [11],  and  Bangs  and  Schultheiss 
[12]  to  name  a  few.  The  purpose  of  this  note  is  to  point  out  that  the 
theory  of  multichannel  maximum  likelihood  processor  yields  a  processor 
structure  that  appears  to  be  attractive  for  adaptive  implementation. 

We  shall  show  that  the  structure  shown  in  Figure  2  follows  naturally 
from  the  theory  of  multichannel  maximum  likelihood  processors.  This 
structure  has  the  advantage  that  the  adaptive  part  of  the  processor  is 
determined  by  the  spatial  and  temporal  properties  of  the  noise  field 
and  it  is  common  to  all  beamsteering  directions.  The  adaptive  processor 


is  followed  by  a  beamsteer ing  operator  which  can  simultaneously  form  a 
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Figure  2.  Multibeara  adaptive  space-time  processor 
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set  of  beams  in  all  the  desired  beamsteering  directions.  The  beam¬ 
steering  operator  consists  of  a  set  of  time  delays. 

In  this  paper,  lower  case  letters  denote  scalars,  underlined  lower 
case  letters  denote  vectors,  and  underlined  capital  letters  denote 
matrices.  The  superscripts  T,  H,  and  *  denote  transpose,  Hermitian 
transpose,  and  complex  conjugate.  All  integrations  are  carried  out 
over  the  signal  duration  T  of  the  active  reception  interval.  I_  denotes 
identity  matrix  and  6(t-z)  denotes  a  delta  function.  Other  symbols 
are  defined  in  the  text. 

II.  MULTICHANNEL  LIKELIHOOD  RATIO  PROCESSOR 

The  structure  of  a  multibeam  processor  for  active  detection  of 
slowly  fluctuating  point  target  in  anisotropic,  nonstationary,  colored 
Gaussian  noise  is  a  generalization  of  multi-input,  single  output  like¬ 
lihood  ratio  processor.  The  likelihood  function  for  an  input  vector 
r(t)  is 

■  r 

*>1  *  vj(u)  Q(u,t)  r(t)  dt  du 

V  V 

T  T 

where  v^(u)  is  the  desired  signal  vector  arriving  from  the  direction  8^, 

Q  (t,u)  is  the  matrix  inverse  kernel  that  satisfies  the  integral  equation 
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Rn(t,u)  is  the  noise  covariance -matrix  [13][14].  The  vector  of  max¬ 
imum  likelihood  functions  Z_  is  given  by 

A-  v  (u)  Q(t,u)  r (t)  dt j  du  (3) 

_  „ 

T  T 

T 

where  Z_  =  Z^  ...  and  the  matrix  V (u)  is  given  by: 

V(u)  =  [v1(u),  v2(u)  ...  vL(u)]. 

The  inner  integral  in  Equation  3  is  an  optimum  spatial  and  temporal 
smoothing  filter  where  the  kernel  Q(t , u)  is  determined  through  Equation 
1  from  the  noise  covariance  matrix  Rn(t,u).  The  covariance  matrix 
Rn(t,u)  describes  the  spatial  and  temporal  properties  of  the  anisotropic, 
colored,  and  nonstationary  noise  field.  The  very  essence  of  the  adaptive 
processing  is  to  "learn"  the  covariance  matrix  R^Ct.u)  from  the  real 
data  and  to  find  a  suitable  approximation  to  the  inverse  kernel  Q(u,t). 

These  two  steps  are  the  most  difficult  and  computationally  most  demanding 
operations  in  the  adaptive  processing.  Hence,  these  steps  should  be 
done  once  for  all  the  beamsteering  directions.  The  multiple  beamsteering 
and  matched  filtering  operation  is  performed  by  the  operator 

* 

V^(u)  •  du 
T 

A  typical  term  of  the  vector  of  likelihood  function  is 


TV 
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^SL  =  v^(u)-y(uO  du 


=  [vA1  (u)  y1  (u)  +  v£2  (u)  y2  (u) 


...+V£KyK]du 

where  y (u)  is  the  output  vector  of  the  adaptive  matrix  filter.  Equation 
4  indicates  that  a  typical  likelihood  function  is  computed  by  a  vector 
correlation  receiver.  The  meaning  of  the  processor  structure  becomes 
clearer  if  we  convert  the  vector  correlation  receiver  to  the  vector 
matched  filter  form: 


h£(T  ~  u)  y(u)  du 


where 


h^(u)  =  v^(T  -  u) 


If  we  now  make  the  usual  assumption  that  the  received  signals  imping¬ 
ing  on  the  sensors  or  subarrays  are  identical,  except  for  the  relative 
propagation  delays  from  sensor  to  sensor  ,  the  8,-th  received  signal 
vector  becomes: 


Vn  (u)  =  [s*  (u  -  T,  )  .  .  .  S  p  (U  -  X  )  .  .  .S0  (U  -  T„)  ] 


This  is  not  an  overly  restrictuve  assumption  because  one 
can  easily  compensate  for  different  sensor  characteristics 
in  the  beamformer.  Many  problems  of  practical  interest 
fall  into  this  category. 
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or  the  vector  matched  filter  can  be  implemented  by  a  set  of  delays 
that  will  bring  the  £-th  look-direction  signal  into  time  coincidence 
at  the  £-th  output  of  the  beamsteering  operator.  That  is,  the  beam¬ 
steering  operator  forms  a  set  of  uniformly  shaded  beams  in  each  of  the 
desired  look-directions.  As  it  is  shown  in  Figure  2,  the  outputs  of 
the  beamsteering  operator  are  processed  by  a  set  of  scalar  matched 
filters. 

III.  ADAPTIVE  IMPLEMENTATION  OF  THE  LIKELIHOOD  RATIO  PROCESSOR 

From  the  previous  discussion,  it  is  obvious  that  the  matrix  inverse 
kernel  Q(t,u)  is  the  only  part  of  the  likelihood  ratio  (LR)  processor 
that  is  based  on  information  which  is  not  available  a  priori:  the 
matched  filters  are  uniquely  determined  by  the  (known)  signals,  and 
the  beamsteering  is  only  dependent  on  the  desired  look-directions. 

The  matrix  inverse  kernel  Q(t ,u) ,  on  the  other  hand,  is  exclusively 
dependent  on  the  physical  properties  of  the  interference  and  these 
cannot,  in  general,  be  assumed  to  be  known  in  advance.  This  clearly 
suggests  an  adaptive  implementation  of  this  part  of  the  processor. 

In  order  to  derive  an  adaptive  realization  of  the  operator  * ] 
as  defined  by: 

r 

(£r)(t)  =  £(t,u)  .r(u)  du 

T 

where  Q(t,u)  is  the  solution  to  the  integral  Equation  2,  it  is  more 
convenient  to  consider  the  "dual"  problem,  i.e.,  the  operator: 
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(Hr)  (t)  =  [(I-Q)r](t)  £  -  [l6(t-u)-Q(t,u)]r(u)da 

T 

*  H(t  ,u)  _r  (u)  du  (9) 

- 

T 

where  H(t,u)  satisfies  the  integral  equation  (substituting  9  into  2) 

- 

H(t, lOR^u, z)  du  =  Rc(t,z)  (10) 

‘ 

T 

jif 

where  R^t.z)  is  the  correlation  function  of  a  "colored"  noise  component. 

It  can  be  shown  from  Equation  10  that  the  linear  system  H[*]  corre¬ 
sponds  to  a  minimum  mean-square  error  estimate  of  this  "colored"  component 
based  on  observation  of  the  input  signal  £(t)  over  the  observation  inter¬ 
val  (in  ether  words,  so-called  optimal  "smoothing") ,  since  the  "ortho¬ 
gonality"  condition  (Papoulis  [15]): 

E{[c(t)  -  (Hr)  (t)]  •  rH(u)}  =  0  t,u  £  T  (11) 

reduces  to  condition  (10) . 

Various  interpretations  can  now  be  made.  While  the  operator  H[*] 
can  be  implemented  as  a  MMSE  estimator  ("smoother")  using  a  multi- 

* 

It  can  be  shown  (Van  Trees  [13])  that  the  relation  between 
this  "generalized"  colored  component  and  the  usual  defi¬ 
nition  of  the  colored  noise  c'(t)  (where  ^.'(t.u)  does  not 
contain  any  singularities)  is  given  by: 

R^t.u)  =  R^,  (t,u)  +  (N  -  I_)  S(t  -  u) 

where  N  is  the  matrix  of  (cross)  spectral  heights  of  the 
temporally  uncorrelated  noise  components. 
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dimensional  version  of  Widrow1 s  -unconstrained  LMS  algorithm  (Widrow 
[16]),  the  actual  "inverse"  filter  Q[*]  which  is  given  by: 

[fir]  (t)  =  r(t)  -  [Hr]  (t)  =  r(t)  -  c(t)  (12) 

(see  Figure  3)  is  obviously  a  "noise"  canceller  (Widrow,  et  al.  [4]). 

It  is  to  be  noted  that  under  general  interference  conditions,  the 
LMS  algorithm  may  have  unacceptable  convergence  properties  due  to  the 
eigenvalue  spread  of  the  input  correlation  matrix  (Widrow,  et  al.  [17]). 

Other  structures  that  have  been  proposed  more  recently  for  noise  cancelling 
applications  (Griffiths  [18] [19])  are  the  so-called  "lattice"  or  "ladder" 
configurations,  which  seem  to  have  efficient  convergence  behavior  as 
well  as  interesting  roundoff-noise  characteristics  (Crochiere  [20]). 

While  these  structures  have  been  primarily  used  for  linear  prediction, 
the  smoothing  case  is  only  a  straightforward  extension  (Friedlander  [21]; 

Lee,  et  al.  [22]). 

As  an  additional  comment,  we  can  add  that  in  the  adaptive  reali¬ 
zation  of  fi[*],  only  those  data  should  be  used  which  are  known  to  be 
interference  data  only.  Since,  of  course,  it  is  unknown  a  priori  whether 
or  not  a  target  is  present,  those  data  corresponding  to  the  range  in 
question  should  be  excluded  from  the  data  used  in  the  implementation. 

This  obviously  limits  the  application  to  the  case  of  active  sonar  or 
radar. 

Finally,  we  should  mention  that  exactly  the  same  structure  as  in 
Figure  2  can  be  used  for  maximum  likelihood  parameter  estimation  (Van 
Trees  [13]).  Thus,  parameter  estimation  does  not  require  different 
algorithms  nor  a  higher  computational  load. 


aik.y  .A* .  iVi 


44 


V'lTTr 
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IV.  CONCLUSION 

A  structure  for  the  implementation  of  a  likelihood-ratio  detector 
is  derived.  This  structure  consists  of  a  "noise  canceller"  which  can 
be  implemented  by  unconstrained  adaptive  algorithms,  followed  by  a 
fixed  beamformer  and  a  set  of  (scalar)  matched  filters.  This  structure 
seems  to  be  more  interesting  than  previous  applications  where  adaptive 
constrained  algorithms  were  used  as  adaptive  beamformers.  Since  no 
constrained  algorithms  have  to  be  used  and  the  space-time  adaptive 
processor  is  common  to  all  "look  directions",  it  is  suggested  that  this 
implementation  will  be  more  computationally  efficient.  Furthermore, 
since  noise  cancelling  has  been  extensively  studied,  various  forms  for 
the  implementation  of  the  adaptive  part  of  this  likelihood-ratio 
processor  are  available. 
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